Topological Anderson Insulator in Three Dimensions 

H.-M. Guo, 1 * 2 G. Rosenberg, 1 G. Refael, 3 and M. Franz 1 

department of Physics and Astronomy, University of British Columbia, Vancouver, BC, Canada V6T 1Z1 
2 Department of Physics, Capital Normal University, Beijing 100048, China 
'Department of Physics, California Institute of Technology, Pasadena, CA 91125 

(Dated: June 15, 2010) 



When the spin-orbit coupling generates a band 
inversion in a narrow-bandgap semiconductor 
such as Sb^Bii-a; or Bi2Se3 the resulting system 
becomes a strong topological insulator (STI)PU2I. 
A key defining property of a STI are its topologi- 
cally protected metallic surface states. These are 
immune to the effects of non-magnetic disorder 
and form a basis for numerous theoretically pre- 
dicted exotic phenomenaSHZl as well as proposed 
practical applications 89 . Disorder, ubiquitously 
present in solids, is normally detrimental to the 
stability of ordered states of matter. In this let- 
ter we demonstrate that not only is STI robust 
to disorder but, remarkably, under certain condi- 
tions disorder can become fundamentally respon- 
sible for its existence. We show that disorder, 
when sufficiently strong, can transform an ordi- 
nary metal with strong spin-orbit coupling into 
a strong topological 'Anderson' insulator, a new 
topological phase of quantum matter in three di- 
mensions. 

Disorder is well known to play a fundamental role in 
low-dimensional electronic systems, leading to electron 
localization and consequent insulating behavior in the 
time-reversal invariant systems^. Disorder also under- 
lies much of the phenomenology of the integer quantum 
Hall effedPl Recently, in a remarkable development, 
it has been noted first by numerical simulationiPZl and 
shortly thereafter by analytical studies^, that a phase 
similar to the two dimensional topological i nsula tor (also 
known as the quantum spin-Hall insulatoi l 14 * 15 |) can be 
brought about by introducing non-magnetic disorder into 
a 2D metal with strong spin orbit coupling. This new 2D 
topological phase, referred to as topological Anderson in- 
sulator (TAI) , has a disordered insulating bulk with topo- 
logically protected gapless edge states that give rise to 
precisely quantized conductance e 2 jh per edge. In TAI, 
remarkably, conductance quantization owes its very exis- 
tence to disorder. 

A question naturally arises whether such behavior can 
occur in three spatial dimensions. More precisely, one 
may inquire whether an inherently 3D topological phase 
analogous to the strong topological insulator^"^ (STI) 
could be reached by disordering a clean system that is ini- 
tially in a topologically trivial phase. This is a nontrivial 
question because just as the 3D STI cannot be reduced 
to the set of 2D topological insulators, the existence of 
a 3D 'strong' TAI presumably cannot be deduced from 
the physics of 2D TAI. Below, we show the answer to the 



above question to be affirmative. Employing a combina- 
tion of analytical and numerical methods we construct an 
explicit example of a disorder-induced topological phase 
in three dimensions with physical properties analogous 
to those of the strong topological insulator. We propose 
to call this new phase a 'strong topological Anderson in- 
sulator' (STAI). We argue that some of the topologically 
trivial compounds with strong spin-orbit coupling dis- 
cussed in the recent literature, such as e.g. Sb2Se3, could 
become STAI upon introducing disorder. In other com- 
pounds that already are STIs in their clean form, disorder 
can reinforce this behavior by rendering the bulk truly 
insulating. We note that the authors of Ref. [T3] antici- 
pated the existence of a 3D disorder-induced topological 
phase. 

To study the emergence of the STAI we consider a 
variant of a model describing itinerant electrons with 
spin-orbit coupling on a cu bic latt ice discussed exten- 
sively in the recent literatur e! 4 ! 19 * 20 !. It has four electron 
states per lattice site rj, compactly denoted as = 
(ipij,ip2j)^3j,ipij) j and a momentum space Hamilto- 
nian Ho = ^ k ^^riu^u with V^k the Fourier transform 
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H k = ^^(k)r M + d 4 (k)i, (i) 

and do(k) = e — 2t^ i cos/c i , dj(k) = — 2Asinfci (i = 
1, 2, 3) and d 4 (k) = 27^(1 - cosfe*)- Here 1^ are 4 x 4 
Dirac matrices in combined orbital and spin space, satis- 
fying the canonical anticommutation relation {r^, r„} = 
2<5 M „. The system defined by Hq is invariant under 
time-reversal and spatial inversion. In the following we 
take two alternate points of view regarding Hamilto- 
nian ([T]): (i) we view it as a simple toy model conve- 
niently describing both topological and ordinary phases 
of non-interacting electrons in 3D and, (ii) we regard 
it as a lattice regularization of the effective low-energy 
Hamiltonian describing the physics of insulators in the 
Bi 2 Se 3 family^JM] j n tne [ a ^ er interpretation la- 
bels the manifold of four active orbitals (|P1 + , f) , \P2~ , f 
), |Pl+,4), |P2~,4}). In addition, parameters t, A and 7 
in these materials show uniaxial anisotropy. We shall ig- 
nore this feature for now but return to it in our discussion 
of the physical realization of STAI. 

The energy spectrum of Hq has two doubly degenerate 



2 



bands, 



£ k = d 4 (k)± Ew- 
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At half filling, depending on the values of the parameters 
e, t, A, and 7 the system can be a metal, a trivial insula- 
tor, as well as a strong and weak topological insulatoif^fll. 
Below, we focus on the topological phase transition be- 
tween t he or dinary insulator characterized by the Z 2 
invarianpE3 (0;000) and the (1;000) STI phase. In the 
clean system modeled by Hamiltonian ([I]) this transition 
occurs due to the band inversion at the T-point of the 
Brillouin zone and can be driven by tuning parameter e 
through a critical value e c = 6t. We take t, A and 7 pos- 
itive here and in what follows. We note that it is exactly 
this physics that underlies the STI behavior in the Bi2Se3 
family of materials^^I. Our main result is the finding 
that a similar transition can be effected by introducing 
non-magnetic disorder into a system that is topologically 
trivial in its clean form. 

To simulate the effects of disorder we consider a Hamil- 
tonian of the form 



(3) 



where Hq is the disorder-free Hamiltonian discussed 
above and Uj is a random on-site potential uniformly 
distributed in the range (— Uo/2, Uo/2). Following the 
discussion in Ref. [13] we start by treating the disorder 
within the self-consistent Born approximation (SCBA). 
The disorder-averaged electron propagator g(u), k) = 
(uj + iS + Ep — Hk — £k) _1 is then given in terms of disor- 
der self-energy £k, subject to the self-consistent equation 
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Here Ep refers to the Fermi energy and S is a posi- 
tive infinitesimal. The factor of 12 arises from the vari- 
ance (U 2 ) — U§/12. General symmetry consideration re- 
strict the form of the self-energy to £ k = Yin ^^^(k) + 
lS4(k). In addition, for point-like disorder £ k is mo- 
mentum independent. This forces £j(k) to vanish for 
i = 1,2,3 as any non-zero value would signal sponta- 
neous time-reversal symmetry breaking caused by non- 
magnetic impurities. 

The two non-vanishing components of the self-energy, 
£0 and £4, can be viewed as disorder-induced renormal- 
izations of the 'topological mass' m = do(k = 0) = e — 6i 
and the Fermi energy Ep, respectively. Specifically, 
it is easy to see that within the SCBA the disorder- 
averaged system is described by the same propagator 
go(ui, k) = {uj + i5 + Ep — "Hk) _1 as the clean system but 
with parameters to and Ep replaced as to — > fh = m+£o 
and Ep —> Ep = Ep — £4. We note that £ and £4 are 
generally complex- valued, giving m and Ep both real 



and imaginary parts. The latter reflect the quasiparticle 
lifetime broadening due to disorder. 

These considerations underlie the physical picture be- 
hind the emergence of STAI. Starting from a clean or- 
dinary insulator (characterized by m > 0) disorder can 
induce a band inversion by driving the real part of the 
renormalized mass to negative. According to the stan- 
dard classification of topological insulators^ a band in- 
version at an odd number of time-reversal invariant mo- 
menta changes the Zi class of the material. If the renor- 
malized Fermi energy Ep lies within the gap the resulting 
effective medium is a STI. 

It is easy to obtain the self-consistent equations for to 
and Ep from Eq. (El). These read 
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We shall present 
the full solution to these below. To gain some insight 
into the underlying physics it is useful to first study the 
approximate solution valid at weak disorder (small Uq) 
obtained by replacing to and Ep on the right hand side 
by their respective bare values. One then obtains 
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where we kept only the leading divergent terms after ex- 
panding the integrand around the T point. We observe 
that for t > 7 disorder indeed renormalizes the mass term 
downward; for Uq > U c ~ [24-7rm(t 2 - 7 2 )/i] 1 / 2 band in- 
version occurs and one expects the system to become 
a STAI. We demonstrate below that this conclusion re- 
mains valid when Eqs. (5j6| are solved self-consistcntly, 
as well as in the numerical simulation of Hamiltonian (|3| . 

The expression for the critical disorder strength U c de- 
rived above suggests that a significant amount of disor- 
der might be necessary to effect the band inversion in a 
real material. Since the Born approximation is expected 
to hold in the limit of weak disorder it is important to 
establish the stability of STAI phase by means of a com- 
plementary method. To this end we carry out exact nu- 
merical diagonalization studies of our model Hamiltonian 
([3]). Generalizing the approach of Refs. [12113125] we 
study the emergence of topologically protected gapless 
surface states that are the defining feature of the strong 
topological insulator in 3D. In a disordered system it is 
not sufficient to establish the existence of gapless modes 
as these could be localized in space. Instead, one must 
seek extended states capable of carrying electrical cur- 
rent across the macroscopic sample. With this in mind 
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FIG. 1: Conductance in disorder-induced phases, a) Conductance G of a rectangular wire with dimensions 4 x 20 x 40 as 
a function of disorder strength Uq for periodic boundary conditions along x and periodic (open) boundary conditions along y as 
distinguished by blue (red) symbols, with Ef = 20meV. Error bars reflect the conductance fluctuations in the ensemble of 100 
independent disorder realizations, b) Electrical current distribution in the STAI wire averaged over 100 independent disorder 
realization for Ef = 20meV and Uq = 150meV. We use open boundary conditions along y and periodic along x direction. 
Arrows representing the local current density in different layers along the y-direction have been color-coded for clarity, c) False 
color plot of conductance as a function of disorder strength Uo and the Fermi energy Ef- Each data point corresponds to a 
single realization of the disorder potential. The color scale is chosen to emphasize the effect of fluctuations in G. The dashed 
line marks the band inversion boundary defined as Re(m) = based on the SCBA. Solid lines represent the SCBA phase 
boundaries separating a band insulator and a metal defined by |Re(-Ep)| = [Re(m)|. d) Density of states at the Fermi energy 
D(Ef) calculated using SCBA. In all four panels we use parameters e = 145meV, t = 24meV, A = 20meV, and 7 = 16meV, 
corresponding to m = lmeV. 



we employ the recursive Green's function method^ to 
evaluate the conductance G of a length-!/ wire with a 
rectangular cross- sectio n W x x W y using the Landauer- 
Biittiker formalis m 1 27 ! 28 !. For simplicity and concreteness, 
and to make contact with the previous works on 2D 
systems ^ 2 * 13 ! 25 ! , we use model parameters listed in the 
caption of Fig. [I] with va lues c lose to those describing 
HgTe/CdTe quantum wellsPUSI 

Figure [T^, shows conductance G as a function of disor- 
der strength Uq in a wire with W x = 4 and W y = 20. For 
weak disorder the wire shows conductance characteris- 
tic of a disordered metal with significant fluctuations re- 
flecting different realizations of the disorder potential £/j. 
Above Uq ~ 130meV the bulk conductance (measured 
with periodic boundary conditions along x and y) drops 
to zero, indicating a disorder-induced insulating behav- 



ior in the bulk that persists up to Uq ~ 170meV. If we 
change periodic boundary conditions to open along the 
y-direction a very different picture emerges. For Uq in 
the range showing bulk insulating behavior the conduc- 
tance is now pinned to the non-zero value 2e 2 /h with no 
observable fluctuations. We attribute this to the ballistic 
transport of topologically protected surface states. Fig- 
ure [TJd confirms that the current flows near the surface 
only. It furthermore shows that the current density is 
equally distributed among the two surfaces; each surface 
forms an independent conduction channel contributing 
one quantum e 2 /h to the total conductance G of the 
sample. This is the disorder-induced topological phase 
mentioned above. 

How can one be sure that this is a genuinely 3D topo- 
logical phase analogous to STI and not merely a result of 
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2D physics already discussed in Refs. [12113] ? First, we 
note that our model Hamiltonian is completely isotropic 
in 3D space so it is unlikely that the observed behav- 
ior would originate from the edge states of a set of 2D 
layers. We checked that the conductance remains quan- 
tized at 2e 2 /h when we change the number of layers W x 
and there is no even/odd effect that one would expect in 
a layered system. The odd number of ID conductance 
channels per surface reflects the odd number of gapless 
states per surface of a STAI, confined to a ID geometry 
imposed by the finite width of the wire in the x-direction. 
When we impose anti-periodic boundary conditions for 
the electron wavefunctions along the x-direction the con- 
ductance drops to zero in accord with the expectation for 
the surface state in a 3D topological phase. 

To further confirm the 3D nature of the observed topo- 
logical phase we probed for the Witten effeclP^ in our 
model STAI. According to Refs. |4|5] the effective electro- 
magnetic Lagrangian of a 3D strong topological insulator 
contains an unusual 'axion' term ~ 6*E • B with 9 = tt. 
According to WitterP^ a magnetic monopole inserted 
into a medium with non-zero 9 binds electric charge 
—e(9/2ir + n) with n integer. Using numerical methods 
described in Ref. [2U] we measured the induced fractional 
charge in a configuration containing a monopole and an 
anti-monopole depicted in Fig. [2^,. Our results presented 
in Fig.[2}3,c clearly indicate fractional charge ±e/2 bound 
to the monopole, confirming the expected value of 9 = n. 
This result lends additional support to our identification 
of STAI as a genuinely 3D topological phase character- 
ized by the bulk axion term. 

We now turn to the phase diagram of our model. 
Specifically, we wish to map out the locus of points in the 
space of parameters (Ep, U ) that gives rise to STAI be- 
havior. In a real material these parameters can be tuned, 
at least in principle, by adjusting the chemical composi- 
tion and disorder content. Because of the 3D nature of 
our system and the resulting large size of the Hamilto- 
nian matrix that must be diagonalized, we were able to 
consider only a single realization of the disorder potential 
Ui at each point of the (Ep, U ) phase diagram. Never- 
theless this turns out to be sufficient for determining the 
location of STAI phase to a good accuracy. Our method 
relies on the fact that, as seen in Fig. [T^,, the conduc- 
tance G shows no observable fluctuations in the STAI 
phase but fluctuates significantly elsewhere. Figure [TJ; 
displays G in a fashion that is designed to amplify the 
effect of fluctuations. The locus of STAI phase is clearly 
visible as the region with G = 2e 2 /h and no discernible 
conductance fluctuations. 

The weak-disorder boundary of the numerically de- 
termined STAI phase coincides with the metal-insulator 
transition deduced from the SCBA Eqs. (7j8|. As in 2D, 
the weak-disorder STAI phase boundary marks the cross- 
ing of a band edge rather than a mobility edge^ ; . The 
strong-disorder phase boundary is more interesting. Ac- 
cording to SCBA (Fig. [Ti) states appear at the Fermi 
level beyond certain disorder strength U* (Ep) due to life- 
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FIG. 2: Witten effect in the strong topological Ander- 
son insulator. a) A prism-shaped sample with monopole 
(+), anti-monopole (— ), and field lines indicated. The asym- 
metric field distribution allows us to use periodic boundary 
conditions in all directions and thus avoid difficulties asso- 
ciated with the surface states, b) Electric charge density 
5p(z) = ~}2 X y [p(r) — po(r)] induced by the monopole/anti- 
monopole pair in a 5 x 5 x 100 sample at half filling. Here p 
and po represent the charge density with and without the pair, 
respectively. A smoothed charge density function Sp(z), ob- 
tained by convolving Sp(z) with a Gaussian of width a — 1.5, 
is also plotted to emphasize charge bound to the monopole. c) 
Integrated charge SQ(z) = ^2 Z , <Z Sp(z') in units of e for Ep 
in the range 19 — 25meV, corresponding to the STAI phase. 
The number of filled electron states is indicated in the leg- 
end (5000 represents the half-filling) . The steps of magnitude 
±e/2 at the monopole/anti-monopole locations show the ex- 
pected localized fractional charge due to the Witten effect. 
Steps with magnitude ±e located elsewhere are due to the 
shift of some bound states in applied magnetic field and can 
be viewed as a finite-size effect which should diminish for a 
system with a larger cross section. In both panels we use 
Uo = 150meV and other parameters as in Fig. 1. 



time broadening. The band edge becomes ill-defined here 
as the self energy acquires a large imaginary part . Com- 
parison to Fig. [Tj: shows that the numerically determined 
STAI phase extends well beyond U*(Ep). The bulk elec- 
tron states near the Fermi level do not contribute to con- 
duction in this regime and must therefore be localized. It 
is tempting to identify this region as the 'true' topologi- 
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cal Anderson insulator where electron localization plays 
the key role. 

Although discussed here in the framework of a con- 
crete model we expect the emergence of disorder-induced 
topological phases to be quite generic in three spatial 
dimensions. We have verified by an explicit calcula- 
tion that including uniaxial anisotropy characteristic of 
the materials in the Bi2Se3 familjEIHS] leads to qual- 
itatively similar behavior to that displayed in Fig. [T] 
This suggests that clean Sb2Se3, predicted to be a triv- 
ial insulator^ (but nevertheless close to the topologi- 
cal phase), could become STAI upon introducing non- 
magnetic disorder. Whether or not disorder can effect 
a band inversion in Sb2Se3 will depend crucially on the 
magnitude of its native bandgap. Our simulations sug- 
gest that when the bandgap size becomes larger than 
~ 20meV the amount of disorder required to produce 



a band inversion is so large that the electron states be- 
come localized before STAI phase can be reached. From 
this point of view the best prospects for experimental 
realization of STAI physics lie with materials that are 
very-small or zero bandgap semiconductors with strong 
spin-orbit coupling. Although bulk HgTeas well as the 
recently discussed Heusler compound d 30 * 31 ! exhibit this 
type of behavior their low-energy physics is not well de- 
scribed by our four-band model and further theoretical 
work is needed to determine whether disorder could drive 
the transition into the topological phase. 
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